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function z= pshiftl(z,th); 
iz=real(z); qz=imag(z); 

tha=2*(bitand(th+32,2. A (0:7))>0)-l; %make array of phase bits + 32 

%0,90, 180,270 

if tha(7)= =l;z=j*z;end; 

if tha(8)==l;z= -z; end; 

% + - 33.69,11.31 
if tha(6)*tha(5)= = -l; 

z=3/16*(5+tha(6)*j)*z; 

else; 

z=17/64*(3+tha(6)*2*j)*z; 

end; 

% +- 8.53,2.84 

if tha(4)*tha(3)= =-1; 

z=(20*( 1+1/1 28) + tha(4)*j)*z/32; 

else; 

z=(20+tha(4)*3*j)*z/32; 

end; 

% -2.86 -1.43 0 1.43 
if tha(l)==l; 

z=(40+tha(2)*j)*z/32; 
elseif tha(2)= = -l; 

z=(20-j)*z/16; 

else; 

z=5*z/4; 

end; 
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Anatomy of the FFT 



This is a 256 point base 2 example of an FFT (a fast dft) 
k:= 0. 255 FFT index 

^ := rnd( l ) + j md( 1 ) FFT dummy arguement 

The discfrete forrier transform takes 256 A 2 operations 
DFT. " ^ 256 



i = 0 

if we combine terms from the first and 2nd half of the summation we have 

127 / S -2, 



I X. -f- X. 



+ 128 

i = 0 

Note the term in paremphases runs over half of N and is only unique for k mod 2 

SO" 

il :=0.. 127 kl:=0.. 1 
men 127 ■ 2z _i± 

dft v= E X1 .,- (k , 2 ) e 256 Z(I DFT *- DFT1 *I) 2=0 

i = 0 k 

This operation only takes half the DFT steps by taking 128 steps to precompute X1 
We can further reduce the computational load by doing the same thing again 

il :=0.. 63 kl:=0.. 3 



i = 0 k 

FIG. 23a 



We repeat this operation a total of 8 times until i ranges only over zero 
We show one more trick. Instead of: 

il :=0.. 31 kl :=0.. 7 



Use 

kl :=0.. 3 



It only requires half the phase shifts. This operation is called a butterfly 
Continuing 

il :=0.. 15 



and finally 



£(l DF v 

k 

FIG. 23b 
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FIG. 27b 
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% finds likelihood of 1 6ary symbol error for nbe erasures 
tr={ 1 -1 ; 1 1]; tr=[tr -tr; tr tr]; tr=[tr -tr; tr tr]; tr=[tr -tr; tr tr]; 
v=zeros(l,16); 
nsym=zeros( 1,16); 
while sum(v)<16; 

[a,n]=max (tr*v'); 
if n<16; 

nbe=sum(v= =0); 
nsym(nbe)=nsym(nbe) +1; 

end; 
k=l; 

v(k)= -(v(k)-l); % -1 or 0 - erasures 
while v(k)= =0; 
k=k+l; 

v(k)= -(v(k)-l); % -1 or 0 - erasures 

end; 

end; 

pe=005:.005:.5; 
for k=l:100; 

we(k)=sum(pe(k). A (l : 16) *( l-pe(k)). A (l 5:-l :0).*nsym); 

end; 

rhe=l-((l-pe). A 4+4*pe.*(l-pe). A 3) A 4; % repeat sym hard error 
rea=l-(l-pe. A 4).64; % repeat sym erasure error 
plot (pe,we,pe,rea) 



FIG. 28 
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FIG. 35f-3 
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